Nonequilibrium Spin-transfer Torque in SFNFS Junctions 
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Abstract. We report theoretical results for the nonequilibrium spin current and spin-transfer torque in voltage-biased SFNFS 
Josephson structures. The subharmonic gap structures and high voltage asymptotic behaviors of the dc and ac components of 
the spin current are analyzed and related to the spin-dependent inelastic scattering of quasiparticles at both F layers. 
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FIGURE 1. Schematic of SFNFS junction. The magnetiza- 
tion p, a is along the z axis (the quantization axis for spin), while 
/I/, is at polar angle \jf in the yz plane. The x axis is along 
Ji a x pf,. Notice our choice of axis differs from Ref. [6]. 



Spin-polarized current passing through a ferromagnet 
can transfer spin angular momentum to the ferromag- 
net and exert a torque on the magnetic moment [1, 2]. 
This mechanism offers unique opportunities to manipu- 
late the magnetic state of nanomagnets. Experimentally 
spin-transfer torque driven magnetization precession and 
magnetization reversal have been observed in magnetic 
multilayers [3, 4]. A well studied multilayer system is 
the magnetic nanopillar which consists of a ferromagnet- 
normal metal-ferromagnet (FNF) trilayer connected to 
normal metal electrodes. The typical thickness of each 
layer is several nm, and the diameter of the pillar is of 
the order of 50 nm [4]. When sandwiched between su- 
perconducting electrodes, the FNF trilayer can mediate 
finite Josephson coupling to form a SFNFS Josephson 
junction [5]. In such junctions, scattering of quasiparticle 
at the magnetic interfaces is sensitive to the phase shift in 
each F layer as well as the condensate phase difference 
(j) across the junction. This indicates the spin momentum 
transfer between the quasiparticles and ferromagnets can 
be tuned by varying <j>. Waintal and Brouwer calculated 
the phase-sensitive equilibrium torque in SFNFS junc- 
tions [6]. They also showed the nonequilibrium torque in 
NFNFS junctions acquires novel features, e.g. it can fa- 
vor perpendicular configuration of the two moments [7]. 



Here we investigate the nonequilibrium spin-transfer 
torque in SFNFS point contacts under bias voltage V. 
The setup and the coordinate system are shown in Fig. 
1 . The two F layers are labelled by indices a and b, re- 
spectively. The N spacer is assumed to be transparent and 
in the clean limit. For simplicity we model each F layer 
as a delta function barrier with spin-dependent transmis- 
sion probability D-t ^ Di , so the spin mixing angle [8] 
# = arcsin ^JD~\ — arcsin y^5J. The scattering matrix of 
¥" and ¥ b are related by spin rotation of angle y. Since 
the condensate phase difference evolves at Josephson fre- 
quency C0j — 2eV/h, the spin current in each region, I, 
(i = 1,2,3), also oscillates with time. For example, 

= Iq + £ [l k , c cos(kcojt) +l u sin(kcojt)]. (1) 

k=l 

We use the quasiclassical theory of superconductivity 
[9] and solve numerically the transport equations for 
the Green's functions in each region which obey proper 
boundary conditions at interface F" and F b [8]. The spin 
currents are then determined from the local Keldysh 
Green's functions. Within this approach the proximity 
effect and the multiple Andreev reflection (MAR) at both 
interfaces are fully taken into account. 

The spin current in the N layer, measured in unit of 
NfVfAAh/2 where A is the contact cross sectional area 
and A is the superconducting gap, is shown in Fig. 2 
for \j/ = n/2 and zero temperature. As in the case of 
nanopillars with normal metal electrodes, the dc spin 
current flow is a consequence of spin filtering at the F 
layers. Its voltage dependence, however, is nonlinear and 
possesses subharmonic gap structures due to the onset 
and resonances of MAR processes. At V = 0, the dc 
spin current vector is along x direction, which tends to 
cause the moments to precess around each other [6]. 
As V is increased, the magnitude of dc currents along 
y and z direction, which tend to bring the moments 
towards or away from each other, grow with V. The 
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FIGURE 2. Dc (top panel), cosine (middle panel) and sine 
(bottom panel) part of the first Fourier components of the spin 
current in the N layer. Each panel features the spin current 
along x, y, and z direction. D* = 0.81, Dy = 0.64, iff = n/2. 

ac spin current originates from the interference between 
MAR processes of different order. Its magnitude changes 
rapidly at voltages below 2A and decays to zero at high 
voltages. 

The spin-transfer torque on F" is given by T„(f) = 
Ii (f ) — 12(f)' an d the torque on F b is given by %{t) = 
12(f) —13(f). It is convenient to expand T a /b(t) in Fourier 
series similar to Eq. (1). Fig. 3 shows the k = and k = 1 
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FIGURE 3. Spin-transfer torque on F* in x (upper panel) and 
z (lower panel) direction. Dt = 0.95, D\ = 0.6, \j/ = 7t/2. 

components of %, in unit of NfV/AAh/2, as functions of 
the bias voltage at T = 0. Notice % is perpendicular to jit, 
and lies in the xz plane (\j/—7t /2). An interesting feature 
of the dc torque along x direction is that it changes sign 
around voltage 1.8A/e. The magnitude of the dc torque 
becomes linear to V at high voltages. The ac torque be- 
comes vanishingly small compared to the dc torque for 
V > 2A because the number of subgap Andreev reflec- 
tions, which constitute the dominant contributions to the 
ac spin current, is inversely proportional to V. 
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